A simple approach to approximate quantum error correction based on the transpose 

channel 
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We demonstrate that there exists a universal, near-optimal recovery map — the transpose 
channel — for approximate quantum error-correcting codes, where optimality is defined using the 
worst-case fidelity. Using the transpose channel, we provide an alternative interpretation of the 
standard quantum error correction (QEC) conditions, and generalize them to a set of conditions for 
approximate QEC (AQEC) codes. This forms the basis of a simple algorithm for finding AQEC 
codes. Our analytical approach is a departure from earlier work relying on exhaustive numerical 
search for the optimal recovery map, with optimality defined based on entanglement fidelity. For the 
practically useful case of codes encoding a single qubit of information, our algorithm is particularly 
easy to implement. 



I. INTRODUCTION 

Quantum error correction (QEC) is one of the corner- 
stones of quantum information and quantum computing. 
Since quantum effects are extremely fragile and suscep- 
tible to damage by environmental noise, many quantum 
communication or computational tasks will be impossible 
without the use of QEC to protect the information from 
noise. QEC is thus critical for the success of quantum 
technologies. The idea behind QEC is a simple one — 
information is stored in a particular part of the system 
Hilbert space, cleverly chosen based on the noise process, 
so that a recovery operation can be applied to retrieve 
the information affected by the noise. 

Much of the discussion in the past on error correction 
focuses on perfect QEC, where the recovery operation ei- 
ther perfectly corrects the full CPTP noise channel, or 
perfectly corrects the errors conditioned on the fact that 
fewer than some t errors occurred. However, an example 
of a code designed for correcting errors affected by weak 
amplitude damping noise presented in [l| suggests that 
the requirement for perfect recovery may be too stringent 
for certain tasks. While the smallest known perfect QEC 
code requires at least five qubits to encode a single qubit, 
the code in [l[ uses only four qubits to achieve compara- 
ble fidelity. This illustrates a key advantage of relaxing 
the requirement for perfect QEC — one might be able to 
encode the same amount of information into fewer qubits 
while retaining a nearly identical level of protection from 
the noise process. The four-qubit code is also specially 
designed for the channel in question, a departure from 
standard QEC codes that seek to perfectly correct up to 
t arbitrary errors on the system. This adaptation of the 
code to the noise channel, an idea emphasized later in 
0, is a crucial factor behind the success of the 4-qubit 
code. Such approximate QEC (AQEC) codes reveal the 
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possibility of designing codes that are better tailored to 
the particular information processing task at hand. 

The analysis in [l[ was based on small perturbations 
of the perfect QEC conditions central to the standard 
theory of error correction. Subsequent work on AQEC 
adopted an alternate approach by recasting it as an opti- 
mization problem. One can formulate AQEC as the prob- 
lem of finding the optimal encoding and recovery maps, 
given a noise channel and the information we want to 
encode (qubit or higher-dimensional object), with opti- 
mality defined in terms of a chosen measure of fidelity. In 
this paper, optimality is measured in terms of the worst- 
case fidelity, i.e., the fidelity between the input state and 
the state after noise and recovery, minimized over all pos- 
sible input states, for given encoding and recovery maps. 
This is a triple-optimization problem since one needs to 
optimize over all possible encodings, recovery maps and 
input states. 

The simplest approach to solving this optimization 
problem is to hold either the encoding or the recovery 
map fixed, and then perform the optimization over the 
remaining two variables — the recovery or the encoding 
map, and the input state. The problem can be further 
simplified by looking instead at measures based on entan- 
glement fidelity Q , and characterize the performance of 
the code averaged over some input ensemble. This elim- 
inates the minimization over all input states required for 
the worst-case fidelity measure. The task of finding the 
optimal encoding or recovery map is then numerically 
tractable via convex-optimization methods , but the 
resulting recovery is now optimal for an averaged mea- 
sure of fidelity. Recovery maps which are near-optimal 
for the average entanglement fidelity have also been con- 
structed analytically, first in Q, and more recently in 

For many communication or computational tasks, how- 
ever, one would prefer an assurance that all the informa- 
tion stored in the code is well-protected. In such cases, 
the worst-case fidelity is the appropriate measure for de- 
termining the optimality of encoding and recovery maps. 
The resulting double-optimization problem for a given 
encoding map was examined using semidefinite program- 



2 



ming in (llj . This method however requires a relaxation 
of one of the constraints in the problem, so the recovery 
map found is typically suboptimal. Furthermore, the nu- 
merically computed recovery map is difficult to describe 
and understand analytically. 

In this paper, using the worst-case fidelity measure to 
define optimality, and assuming a fixed encoding, we con- 
struct a universal recovery map that is very easy to define 
analytically. This universal recovery map — the transpose 
channel 0, LL2| — gives a worst-case fidelity that cannot 
be too far from that of the optimal recovery. Using 
the fact that the transpose channel is the optimal re- 
covery map for perfect QEC codes, we rewrite the error 
correction conditions [llj| - [l5j for perfect QEC in such a 
way that the role of the transpose channel is apparent. 
From this, we derive necessary and sufficient conditions 
for AQEC founded upon the transpose channel, as a nat- 
ural generalization of the perfect QEC conditions. While 
AQEC conditions have been derived in the past from 
an information-theoretic perspective [f6l - [20j . our condi- 
tions are algebraic, and lead to a simple and universal 
algorithm to find AQEC codes that does not require op- 
timizing over all recovery maps for each encoding map. 
Furthermore, we demonstrate that the worst-case fidelity 
for the transpose channel is an easily computable quan- 
tity for the most practically useful case of codes encoding 
a single qubit. 

Note that AQEC based on the worst-case entangle- 
ment fidelity was also discussed recently in [21], around 
the same time this work was done. 



II. AQEC AS AN OPTIMIZATION PROBLEM 

Consider a physical system, with Hilbert space de- 
noted by Ti. In this system, we seek to encode a qudit 
of information — information carried by a d-dimensional 
Hilbert space Tia, with d < dim('H). In particular, we 
focus on the case of a subspace code, where the qudit 
is encoded into a d-dimensional subspace C, of Ti. For- 
mally, the information is encoded into C via a linear, 
invertible encoding map W. The action of noise on the 
system is described by a completely positive (CP), trace- 
preserving (TP) map £ : B(H) —> B(Ti). £ can describe, 
for example, the Markovian noise acting on the system 
over some time step, or the effects of a single use of a 
noisy channel for communication. Complete positivity 
of £ entails that its action can be described by a (non- 
unique) set of Kraus operators {Ei}f =1 , such that £ acts 
as £(p) = 5^ i=1 EipEj. To denote the noise channel in 
terms of its Kraus elements, we write £ ~ {Ei}. The fact 
that £ is TP is enforced by the condition ^iEjEi = I, 
where I is the identity operator for the domain of £ . Af- 
ter the action of £, we perform a CPTP recovery map 
TZ : B{Ti) B(C) to undo the effects of the noise, and 
then decode using W _1 . 

How well the information is protected from the noise 
can be quantified by the fidelity between the input qu- 



dit state and the decoded state after noise and recovery. 
The fidelity between any two states p and tr is given by 
F{p,a) = tr-y/ p 1 / 2 o-p 1 / 2 . For a pure state p = \ijj)(ip\, 
this can be written as F(\ip),a) = F (\ifj)(ip\, a) = 

a/OAIoM- Note that, < F(p,a) < 1, with F = if 
and only if p and a have orthogonal support, and F = 1 
if and only if p — a. The fidelity is thus a measure of how 
close two states are. Since we will often discuss fidelity 
for a state before and after the action of a map 3>, we use 
the shorthand F(\ip),$) = F[\ip),$(\ip){il}\)]. 

Based on the fidelity measure, we say that a code C, 
together with W and TZ, is effective at protecting the 
information from the noise £ if the worst-case fidelity 
mm p65(-Ho) ° TZ o £ o W] is close to 1. Here, 

S(Ho) denotes the set of all states, pure or mixed, in the 
codespace. In fact, since the fidelity F{p,u) is jointly 
concave in its arguments (see, for example, [22(), it suf- 
fices to minimize over pure states in S(Jio) only. 

Above, we considered a given encoding map W and a 
given recovery map TZ. In reality, one wants to maximize 
the error correction capability provided by the system by 
choosing W and TZ such that the worst-case fidelity is as 
close to 1 as possible. The problem of AQEC using a 
system with Hilbert space Ti can thus be phrased as 

max max min F(\ip), W" 1 o TZ o £ o W). (1) 

W n |v>e« 

If the quantity in Eq. ([T]) attains the maximum possible 
value of 1, i.e., there exist W and TZ such that the worst- 
case fidelity is 1, then we have perfect QEC. 

The simplest approach to solving this optimization 
problem is to do an exhaustive search over all possi- 
ble encodings. This amounts to randomly choosing d- 
dimensional subspaces C C Ti. For each C, we still need 
to optimize over TZ to maximize worst-case fidelity. For 
a given C, the optimization problem can be written as 



max min F(\ib),TZo £) 

n |V)ec 



(2) 



where the worst-case fidelity is computed over all pure 
states in C only. 

Before proceeding further, let us define some terminol- 
ogy. We will often make use of the square of the fidelity, 
which we denote as F 2 {-,-) = [F(-,-)] 2 . Whenever un- 
ambiguous, we will also refer to F 2 as the fidelity. The 
recovery TZ with the largest worst-case fidelity for a given 
C is the optimal recovery and is denoted by TZ op . The fi- 
delity loss rj-ji. for a given code C and a recovery TZ, is 
defined as 



rm = 1 — min F' 

Wee 



[|V>,(fc°£)(MM)] 



(3) 



The fidelity loss for TZ op , denoted r] op , is 77 op = min^ rm . 
We refer to rj op as the optimal fidelity loss. A code C is 
said to be e- correctable if it has 77 op < e for some e S 
[0, 1]. e-correctable codes with e < 1 are approximately 
correctable in the sense that code states have fidelity at 
least \/l — e ~ 1 — e/2 after the action of the noise and 
(optimal) recovery. 
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III. TRANSPOSE CHANNEL AS UNIVERSAL, 
NEAR-OPTIMAL RECOVERY 

Here, we describe the transpose channel and demon- 
strate that it is indeed the standard recovery map for 
perfect QEC codes characterized by the well-known QEC 
conditions. We then proceed to show that the transpose 
channel is nearly optimal even in the case of AQEC codes. 

A. Transpose channel 

Consider a d-dimensional code C, and a CPTP noise 
channel £ ~ {Ei]f =1 . Let P be the projector onto C and 
Ps be the projector onto Ve = the support of £{C). The 
transpose channel is the CPTP map H P : B{P £ ) -> B(C) 
such that 

N 

TIp (•) = Y, PE l£{.P)~ 1/2 {:)£{P)- 1/2 EiP, (4) 

i=l 

where the inverse of £(P) is taken on its support. The 
transpose channel can be understood as being composed 
of three CP maps: Hp = V o £^ o J\f ', where V is the 
projection P(-)P onto C, £^ is the adjoint of £, and Af is 
the normalization map Af(-) = £(Py 1/2 (-)£(Py 1/2 . In 
this form, Hp is manifestly independent of the choice of 
Kraus representation for £. 

Hp is a special case of a recovery map introduced in Q 
for reversing the effects of a quantum channel on a given 
initial state. Hp defined here is exactly the case for the 
initial state P/d, where d is the dimension of C. In [23| . 
Hp was shown to be useful for correcting information 
carried by codes preserved according to an operationally 
motivated notion. The term transpose channel owes its 
origin to (l2l |. where this channel was first defined in an 
information-theoretic context. It was shown [24[ that the 
transpose channel has the property of being the unique 
noise channel that saturates Uhlmann's theorem i.e. the 
monotonicity of relative entropy — a fact that was later 
used to characterize states that saturate the strong sub- 
badditivity of quantum entropy [25j | . 

While our focus is on AQEC, understanding the rele- 
vance of Hp to perfect QEC provides the intuition be- 
hind the AQEC conditions presented later. An important 
characterization of perfect QEC codes is the set of QEC 
conditions J3-[l5j|, which we briefly review here (see, for 
example, [22|| ) : 

Theorem 1 (Perfect QEC conditions). A CPTP re- 
covery H that perfectly corrects a CP map £ on a sub- 
space code C exists if and only if 

Vi,j, PEjE j P = a ij P, (5) 
for some complex matrix a. 

It is useful to rewrite Eq. §5§ in a "diagonal" form, a is 
clearly Hcrmitian, and can be diagonalized with a unitary 



u such that a = uDv) , where D is the diagonal matrix 
of eigenvalues. The set {F k = J2i u ikEi} constitutes a 
different Kraus representation for £. With this choice of 
Kraus representation, the perfect QEC conditions take 
the diagonal form 

VM, PFtFiP = 5 kl d kk P, (6) 

where d kk are the diagonal entries of D. 

The recovery map H when Eq. is satisfied — which 
we denote as H pe ii — is constructed as follows [22|: us- 
ing the polar decomposition F k P — \/d kk U k P, H pcr f : 

B(V £ ) -> B(C) is given by H pcl - { ~ {PUl). One can 
check that 72.p er f is TP on its domain B(Ps), and that 

for any p G B(C), (H pcli o £)(p) = \YL k d k i?jP- From the 

QEC conditions (Eq. ©), we see that J^k dkk = ^[£{p)\ 
is independent of p, and is exactly equal to 1 if and only 
if £ is TP on C. H pe rf thus recovers the original code 
state, up to any reduction in trace due to the possible 
non-TP nature of £ . 

A natural question to ask here is how the transpose 
channel Hp relates to the recovery H pc ^ for a given £ 
and C that satisfy the QEC conditions. Here, we show 
that they are exactly the same map, as previously noted 
in @: 

Lemma 2. Hp = H per f. 

Proof. Observe that £{P) = Ek( F k p )( PF l) = 
Efe d kkPk, where P k = U k PU%. Eq. © gives PU%U,P = 
SkiP, so that Pfe's are orthogonal projectors with PkPi = 
SkiPk- Hence, £{P)~ 1 ' 2 = Y, k p k/Vd^, The Kraus op- 
erators {PF^,£{P)^ 1 / 2 } of Hp can hence be written as 

PFl^P)- 1 ' 2 = J2 y/d kk /d u PUlUiPUl = PUl (7) 

which are exactly the Kraus operators of 7?. por f. I 

Perfect QEC is often discussed for a noise channel that 
is CP but not necessarily TP. In fact, Theorem [T] and 
Lemma [2] remain true even for a non-TP £. The non- 
TP scenario is particularly relevant when we deal with 
a system of n quantum registers, where each register is 
independently affected by some CPTP noise £\. One 
often looks for codes that perfectly correct the noise up 
to some maximum number t of quantum registers with 
errors. Then, instead of having £ = £f n , the relevant 
noise channel for perfect QEC describes noise where at 
most t registers have errors. Such an £ is not TP, since 
we have discarded the part of £f n that corresponds to 
having errors in more than t registers. 

Actually, a perfect QEC code for such a non-TP noise 
channel can be viewed as an AQEC code for the original 
n-register noise channel £f" 1 , which is TP. In our AQEC 
discussion, the code we look for is approximately cor- 
rectable on the channel anyway, so £ is always assumed 
to be TP, which is often the physically-relevant scenario. 
The TP requirement is also important for fidelity to be 
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a good measure of the efficacy of the recovery operation. 
Note that the analysis in the remainder of the paper ap- 
plies to a special type of non-TP maps — £ ~ {Ei} such 
that Y.i PE \ E i p = O-P for < o < 1, giving an addi- 
tional proportionality factor a in our expressions. 



B. Near-optimality of the transpose channel 

In general, TZp need not be the optimal recovery map 
TZ op for a given C and £ . However, in the following the- 
orem and the subsequent corollary, which form the core 
results of our paper, we show that it does not do much 
worse than lZ op . 

Theorem 3. Consider a d- dimensional code C with opti- 
mal fidelity loss rj op under a CPTP noise channels. For 
any \ip) G C, 



F 2 {\^),K op oE) < ^l + (d-l)r, op F(\iP),n P o£). (8) 

Proof. Let {Rj} be a set of Kraus operators of lZ op : 
B(P £ ) -> B{C). For any \ip) G C, following [|, we have 



F 2 (\^),TZ op o£) 



(9) 



< 



X|(4£(P)-l/2 jBl )|2)(^| (jR ^ (P) l/2 jR t ) |2) ) 



where (■) denotes expectation value with respect to 
Since K op is TP, we have that £\ \(Rj£ (P) 1 / 2 ^])! 2 < 
(J2 j R j S(P)Rt) = ((1Z op oS)(P)). 

Now, choose a basis {\ipi)}f = i for C with \tpi) = 
Let Pi ee (K op o f)(|Vi)<^J) = Eh««I^X^I, fOT 
coefficients satisfying ^ fc aj^ 



From the definition of 7y op , a-^ = (V>i|A;|V'i) > 1 — ??o P - 
This implies a lfc — ^op, which in turn gives 

a kl — Vop Vfc ^ i. Since \tp) = \ipi) by construe- 
tion, we get Ej\(Rj£(P) 1/2 R})\ 2 < ((P op o £)(P)) < 
l + (d— l)??op- Putting this into Eq. ©, and noting that 



(0 
hi 



X.\(E}£(P)- 1/2 E. t 



,1/2 



< F(\i/)),n P o £), gives Eq. 



Let 7yp be the fidelity loss for code C with Rp as the 
recovery. Then, Theorem [3] implies 

Corollary 4. rj P satisfies rj op < n P < i] op f(v op ;d), 
where f(n; d) is the function 

f( m d) EE ( ^ 1} ~ 7? = (d + 1) + 0( V ). (10) 

1 + [a — 1)7] 

Proof. r)p > ?7op is true by definition of n op . For any l^) G 
C, let F 2 (| f),n P o = 1 - r) P ^. Then, by definition, 
the fidelity- loss is r\p = max^(?7p j ^). From Theorem [31 
l-r? op < F 2 {\^),P op o£) < y/ll + id-Vrhplil-ripj). 
Rearranging gives np^ < r\ Qp f{j] op \£). Since this holds 
for all ??p,0, it also holds for r/p. M 



The inequality np < ri op f (r] op ; d) makes precise our 
statement that IZp is near-optimal, with the additional 
factor of (d + 1). For the most practically-relevant case 
of a code encoding a single qubit, this is only a factor of 
3. Note that, for n op = 0, the inequality in Corollary [4] 
collapses to r\p = ry op , as expected from Lemma[5J Corol- 
lary U provides necessary and sufficient conditions for C 
to be approximately correctable — C is approximately cor- 
rectable if and only if rjp is small. 

We do not know if the upper bound on rjp in Corollary 
H]is tight. However, the appearance of the dimension d 
of the code in the bound is unavoidable, as can be seen 
from the following example: 

Example 1. Consider a noise channel £, whose ac- 
tion on a code C is given by the set of Kraus operators 

{EiP} = {VT~^ p,,/p \o){o\,Vp |o)(i|, ...,Vp\o)(d- 

1|} ; for 0<p<l. £ acts like the identity channel on C, 
except for a small damaging component that maps a small 
part of every code state onto |0). For d > 3, one can show 
that the worst- case fidelity, when using IZp as the recov- 
ery, occurs for state |0). The corresponding fidelity loss 
is r/p = (d— l)p/[\ + (d— l)p\. Since £ is nearly the iden- 
tity channel, we might instead do nothing (identity chan- 
nel as the recovery), for which the fidelity loss is r/Q = p. 
rjo is always smaller than rjp for small p. Since r) op < t]q, 
we see that np/i] op > i]p/rj = (d - 1)/[1 + (d — l)p], 
which grows as d increases, for fixed p. Hence, for this 
noise channel and code, the separation between rjp and 
flop grows as d increases. 

That the dimension of the code space appears here is 
perhaps not surprising. In the next section, we will see 
that this approach to AQEC using the transpose channel 
can be viewed as a perturbation from the perfect QEC 
case. The factor of d appearing in our bounds can hence 
be understood as quantifying the number of degrees of 
freedom in which the approximate case can deviate from 
the perfect case. 

Note, however, that as d gets large, f(rj; d) approaches 
1/eta. In this case, the inequality in Corollary 4 simply 
becomes the trivial statement n p < np < 1. While we 
will often only be interested in codes with small values of 
d, this demonstrates the weakness in the bounds derived 
here for large d values. 



IV. THE TRANSPOSE CHANNEL AND QEC 
CONDITIONS 

One of the key tools in perfect QEC are the perfect 
QEC conditions (Theorem [T]). Conditions characterizing 
AQEC codes would likewise be useful. A natural ap- 
proach is to perturb the perfect QEC conditions to allow 
for small deviations. For example, the four-qubit code for 
the amplitude damping channel in [l| was shown to obey 
a set of perturbed QEC conditions. More recently, [26[ 
examined small perturbations of the perfect QEC condi- 
tions for general CPTP channels. However, the analysis 
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in |26| is often complicated, and one wonders if there is 
a simpler approach using the transpose channel. In this 
section, we discuss such a set of AQEC conditions built 
upon Corollary 2J We begin by first writing down an 
alternate, but equivalent set of perfect QEC conditions 
which highlights the role of the transpose channel: 

Theorem 5 (Alternate form of perfect QEC con- 
ditions). A code C satisfies the perfect QEC conditions 
( Theorem Q]) if and only if it satisfies 

Vi,j, PElZ{P)- 1 l 2 E ] P = p i3 P, (11) 
where /? = y/a, for a is defined in Eq.§o\. 

Proof. For a code C that satisfies the perfect QEC condi- 
tions (Theorem [I}, using Eq. and PU ] k UiP = 8 kl P, 
wc have 



PFl£(Py 1 / 2 F l P = S kr 



(12) 



This diagonal form can be rotated to any other Kraus 
representation using a unitary u so that F k = ^ u^Ei 
and a = uDu'. Defining j3 = y/a gives Eq. (fTTj) , thus 
showing that if a code C satisfies the perfect QEC condi- 
tions, it also satisfies Eq. (fTTj) . 

Conversely, suppose we start from the diagonal form 
of Eq. (fTTj) as in Eq. (p"2|) . which can be accom- 
plished by choosing u so that (3 is diagonal with en- 
tries \fd~kk. Then taking the square root of Eq. (fT2|) 
gives £(P)~ 1 ^ 4 F k P = (dkk) 1 ^ 4 VfeP for some unitary 
Vfc, so that F k P = {d kk ) 1/4 £{P) x / A V k P. Putting this 
into Eq. (|12l) gives PVlViP = S k i P. Furthermore, 
£{P) 1/2 = [J2 k (FkP)(PFl)} 1/2 = Z k Vd^V k PVl Di- 
rect computation then gives PF^FiP = 5 k id kk P, which 
is exactly Eq. ©. Applying an appropriate u to rotate 
to the desired Kraus representation gives Eq. (JSJ) - ■ 

Observe that the left-hand side of Eq. (TIT]) is a Kraus 
operator of IZp o £. Thus, the QEC conditions in The- 
orem [5j and equivalently the original conditions stated 
in Theorem [TJ simply express the fact that C is perfectly 
correctable if and only if lZpo£ oc V . The proportionality 

factor is Pij = Eij «.j = Efc dkk- 

We can now obtain a set of conditions for AQEC by 
perturbing this alternate form of the QEC conditions. 



Theorem 6 (AQEC conditions). Consider a CPTP 
channel £ ~ {Ei}, and a d- dimensional code C with pro- 
jector P. Let Aij G B(C) be traceless operators such that 



PEj^Py^EjP = PijP + Ay, 



(13) 



where ftij G C. Then, for e G [0, 1], 3 77 G [0, 1] given by 



77 = max 

Wee 



such that 



[(mlj^m - iMA«Mi s 



(14) 



(i) C is e-correctable if 77 < e; 

(ii) C is e-correctable only if n < ef(e;d), where f is 
the function defined in Eq. (1101) . 



Proof. The left-hand side of Eq. (JT3J) is a Kraus operator 
oiTipo£. This, along with the TP condition for TZp o £ , 
gives the expression on the right-hand side of Eq. ([14")) 
for r\p. Setting 77 = rjp, conditions (i) and (ii) follow from 
Corollary gj ■ 



Eq. (|14p elucidates how the fidelity loss arises from the 
presence of the operators. If A^ — Vi,j, we have 
perfect QEC. 

The AQEC conditions, like the perfect QEC condi- 
tions, provide a way to check if a code is approximately 
correctable, without requiring knowledge of the optimal 
recovery. More precisely, given a maximum tolerable fi- 
delity loss e for some information processing task at hand, 
one can check if a code C is e-correctable as follows. The 
AQEC conditions instruct us to compute rjp, which can 
be done once we know C and the noise channel £. If 
r/p < e, then C is a good code. If however, r\p violates the 
inequality in condition (ii), we know that C is not good 
enough for our purposes. Of course, there is a gap — 
for ijp taking values e < r/p < ef(e;d), we cannot use 
the conditions to determine if C is within our tolerable 
fidelity loss, but this gap is small for small d. We do 
not know if the gap can be shrunk by replacing ?7p with 
the fidelity loss for a different recovery map than the 
transpose channel, but we believe it is unlikely to vanish 
completely. 

For a general C, the fidelity loss r/p may be difficult to 
compute as it requires a maximization over all states in 
the code space. However, there is a quick way to check 
for sufficiency by relaxing condition (i) of Theorem [6j 



Corollary 7. C is e-correctable for 



where A, 



the operator norm. 



Y,,, XA-j- and 



e G [0, 1] if 
denotes 



Proof. Observe that the right-hand side of Eq. (fl"4"l) sat- 

isfics E^IAj.A^) ~ l(V;|Al,#)| 2 ] < MA sum |# 
Maximizing this expression over all G C gives || A sum ||. 
Hence, ijp < ||A sum ||, and the sufficiency condition (i) in 
Theorem[6]is satisfied if ||A sum || < e. ■ 

Since A sum > 0, its operator norm is given by its maxi- 
mum eigenvalue, which is easily computable. In fact, for 
codes encoding a single qubit, it is easy to show (using the 



Pauli basis, for example) that ||A S 



1 - E 47 I/?, 



Note that (3ij for any code C and noise channel £ is simply 
given by (3 t3 = (1/d) tr^^P)- 1 / 2 E 3 P). 

V. COMPUTING 77P FOR QUBIT CODES 

Computing 77^ for a general code requires an exhaus- 
tive optimization over all states in the code. However, 
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for the practically relevant case of codes encoding a sin- 
gle qubit, i.e., C with dimension d = 2, rjp turns out to 
require only simple eigen-analysis to compute. 

For a qubit code, (K P o £ ) : B(C) -> B(C) is a qubit 
map. Observe that IZp o £ is not only CPTP but also 
unital (i.e., (IZp o£)(P) = P). Here, we show that the 
worst-case fidelity for a unital, CPTP qubit map is easy 
to compute. While our context requires only a unital, 
CPTP qubit map, we begin with a general CP map <& ~ 
{Ki} on a d-dimensional Hilbert subspace C, so as to 
highlight why the qubit case is particularly simple. 

We begin by choosing a Hermitian basis 
{O 0> Oi,...,O cP _ 1 } for B(C) where O = I, Ot = 
O q Va, trjO^O/j} = <5 Q /jrf V a, /3. The operators 
{O a ,a = l,...,d 2 — 1} are clearly traceless. Such a 
basis exists for any d — for example, one can use the set 
of standard generators of the SU (d) group, augmented 
with the identity operator. The action of $ can be 
represented as a matrix M , acting on vectors (operators 
in B(C)) in the Hilbert-Schmidt space, with matrix 
elements 



M aP = -tv{O a $(Op)}. 



(15) 



Since $ is CP and O q 's are Hermitian, M* a p = Ma/3, so 
M is a real matrix. 

Now, the density operator corresponding to any pure 
state in C can be expanded in terms of the Hermitian 
basis as 



s-O) = -s-0, 



(16) 



where s is a real (d 2 — l)-element vector, s = (1, s), O = 
(0\, O2, • • • , 0^-1), and O = (I, O). s is not an arbitrary 
vector, but in general has to obey some constraints in 
order for it to correspond to a pure state. 

Using Eqs. (15]) and (|T6]). the fidelity for a state \ip) £ C 
under the map <I> can be written as F 2 [\tp), ^(IV'XV'I)] = 
2 s T M s, where s is s written as a column vector, and the 
superscript T denotes the transpose. We can rewrite the 
expression for the fidelity using the symmetrized version 



oiM: M 



sym 



;(M + M T ). Observe that s T M sym s 



s M s. Therefore, 



f 2 (|v>),$(|V)(V>I) 



'M s 



(17) 



Finding the worst-case fidelity is hence equivalent to the 
following minimization problem for a real, symmetric ma- 
trix M s 



'sym- 

minimize: 
constraint: 



s T M 



sym S, (18a) 

s corresponds to a pure state. (18b) 



For d > 2, the constraint Eq. (|18b[) is difficult to write 
down. Even if we relax the constraint to include mixed 
states, it is not known in general what s corresponding 
to a (positive, trace- 1) density operator looks like. This 
constrained minimization problem is hence not simple for 
a general d. 



For qubits (d — 2) however, the constraint equation 
is simple to write down. We choose the operator ba- 
sis to be the Pauli basis. Given an orthonormal ba- 
sis {\vi), ^2)} for the qubit code space, the Pauli basis 
{(To = l2> a xt Gyi v%\ can be constructed as 

co = |t>i)(ui| + \v-z){vi\ = h, 

<Jx = \vi){v 2 \ + |U2>(«1|) 

a v = -i'(|vi>(«2| - |«2>(«i|), 
and (T z = \v 1 )(v 1 \-\v 2 )(v2\. (19) 

Eq. (|16[) then corresponds to the Bloch sphere repre- 
sentation of a pure state, with the Bloch vector s = 
(s x ,s y ,s z ) satisfying ||s|| = (si + s 2 y + si) 1 / 2 = 1. The 
constraint Eq. (|18b[) becomes 



constraint: s = (l,s), with ||s|| = 1. 



W) 



The constrained minimization problem can then be 
solved using the Lagrange multiplier method. 

For the case of a CPTP qubit map that is also unital, 
the minimization problem can be further simplified. For 
any CPTP, unital $ (arbitrary d), M takes the form 



(1 



M = 





Vo 



... \ 

r 

/ 



(20) 



The first row comes from the fact that $ is TP, and the 
first column from the fact that $ is unital. T is a (d— 1) x 
(d— 1) real matrix. Defining T sy m = ^(T+T^), Eq. (fP7|) 
can be written as F 2 (\ip), <P(\ip)(ij\) = \(l + s T T Bym s). 
This means that we can equivalently minimize s T 7^ ym s 
instead of the original s T Ms ym s in Eq. (|18a|) . For a 
qubit CPTP, unital $ then, the constrained minimization 
problem becomes 



minimize: s 7~ sym s, 



constraint: 



1. 



(21a) 
(21b) 



This simply tells us to minimize the expectation value of 
% y m with respect to all real unit vectors s. 

Since 7^ ym is real and symmetric, it can be diagonalized 
with an orthogonal matrix Q so that % ym = Q t TdQ, 
where Tb is a real, diagonal matrix of eigenvalues of 
T sym . Then s T T sym s = {Qs) t Td(Qs). Q, being or- 
thogonal, preserves the length of the vector it acts on. 
The minimization problem Eq. (f2"Tj) thus corresponds to 
minimizing the expectation value of To over all real unit 
vectors. As Td is real and diagonal, this minimum expec- 
tation value is exactly the smallest eigenvalue of Td (and 
hence of T sy m), attained by the corresponding eigenvec- 
tor normalized to unit length. Therefore, we see that the 
fidelity loss for a CPTP, unital qubit map $ is given by 



?7$ = 1 - min F 2 (|^),$ 



1 



) n (1 ^min); 



(22) 
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where i m i n is the smallest eigenvalue of T syul correspond- 
ing to the map 3>. Setting $ = TZp o £ o V gives rjp. 
Note that, for $ with a Hermitian-closed set [3(| of 
Kraus operators, as is the case for $ = TZp o £ o V ~ 
{PE\£{Py 1 / 2 E\P}, T is symmetric so that % ym = T. 



VI. EXAMPLE: AMPLITUDE DAMPING 
CHANNEL 



" 0.9 



0.8 



0.75 




no error-correction 
■ ■ [4, 1 ] code, Leung recovery 

— [4,1] code, R p recovery 

— [4,1] code, Fletcher recovery 

— [[5,1,3]] code 

— -random 4-qubit code, R p recovery 

— random 4-qubit code. Id recovery 



0.15 0.2 0.25 0.3 0.35 
damping parameter (y) 



FIG. 1: Codes for the amplitude damping channel, for < 
7 < 0.5. 

As an example to illustrate our discussion so far, let us 
look at the noise channel considered in [l[ — the amplitude 
damping channel. The single-qubit amplitude damping 
channel £ ad is the CPTP channel with Kraus operators 



En = 







1 



VT^7 



and Ei = 




(23) 



written in some qubit basis {|0), |1)}. £ad can be thought 
of as describing energy dissipation for a system where |0) 
is the ground state, and |1) is some excited state. 7 is 
then the probability of a transition from the excited state 
to the ground state. In the absence of any encoding or 
recovery, the worst-case fidelity for a single qubit under- 
going £ad fads off as 1 — 7, as 7 increases (see Fig. [TJ 
line labeled "no error correction"). 

A code that uses four physical qubits to protect a single 
qubit of information against amplitude damping noise 
was constructed by Leung et al. [l|. Assuming that the 
noise acts independently on the qubits, the four-qubit 
noise channel is just four copies of ^ADji-e., £alv The 
four-qubit subspace code constructed in |l| is the span of 
the following two states: 

104 = ^(10000) + 11111)), 

and \1 L ) = 4=(|0011> + |1100)). (24) 
v2 



|0l) and \1l) respectively represent the |0) and |1) states 
of the single qubit of information we want to encode in 
the four-qubit Hilbert space. We denote this code as the 
[4,1] code, where the first entry in the brackets corre- 
sponds to the number of qubits in the system, and the 
second entry is the number of qubits of information en- 
coded in the system. It was shown in [l[ that this code 
satisfies the perfect QEC conditions for £®^>, except for 
small corrections of order j 2 , and hence a recovery op- 
eration similar to 7?. per f can be constructed. We refer to 
this recovery map as the Leung recovery. The worst-case 
fidelity for this code and recovery is plotted as a function 
of 7 in Fig. Q] Clearly, the [4,1] code is able to signifi- 
cantly raise the worst-case fidelity for the encoded qubit 
of information, as compared to the no error correction 
case. 

In the same figure, we also plot the worst-case fidelity 
using the transpose channel TZp as the recovery opera- 
tion, instead of the Leung recovery, for the same [4,1] 
code. We see that using the transpose channel as the re- 
covery map gives a higher fidelity than the original Leung 
recovery. 

For comparison, we also look at a recovery map for 
the [4,1] code constructed by Fletcher et al. in This 
recovery map, which we refer to as the Fletcher recov- 
ery, was originally optimized for an averaged measure of 
fidelity. We instead compute the worst-case fidelity for 
this recovery [3l|, also plotted in Fig. Q] For small val- 
ues of 7, the Fletcher recovery gives the best performance 
compared to the other recovery maps, despite being op- 
timized for an averaged measure of fidelity. However, it 
only does marginally better than the transpose channel 
recovery. 

We also compare the performance of the [4,1] code un- 
der these different recovery maps with that of a code 
that satisfies the perfect QEC conditions. The smallest 
code capable of perfectly correcting an arbitrary error 
on any single qubit, requires five qubits. The relevant 
noise channel now is £ ® D . The five-qubit code HH , 
usually referred to as the [[5,1,3]] code [Hj], satisfies the 
perfect QEC conditions for the CP channel comprising 
only the single-qubit (Pauli) errors in Sff). Using the 
corresponding 7^ p0 rf as the recovery for the [[5,1,3]] code, 
we compute the worst-case fidelity for the noise channel 
for different values of 7. As the plot in Fig. [T]shows, 
the [[5,1,3]] code performs better than the [4,1] code with 
Leung recovery, but the [4,1] code uses one qubit less to 
encode the same amount of information. The [4,1] code 
with the transpose channel as recovery has nearly iden- 
tical worst-case fidelity as the [[5,1,3]] code, while the 
one with Fletcher recovery does slightly better for small 
values of 7. 

These observations clearly demonstrate the benefit of 
going beyond codes described by the perfect QEC con- 
ditions. Furthermore, while the [[5,1,3]] code is capable 
of perfectly correcting an arbitrary single-qubit error in 
a system subjected to any noise channel, the comparison 
with the [4,1] code with various recovery maps clearly 



8 



show the gain that one might achieve by adapting the 
codes and recovery to the noise channel in question. 

Lastly, we also compute the worst-case fidelity for ran- 
domly generated four qubit codes, using the transpose 
channel as the recovery map. Computing F 2 for about 
500 randomly selected codes took less than half an hour 
on a typical laptop computer. We plot the worst-case 
fidelity for the best code in Fig. Q] (line marked "random 
4-qubit code, IZp recovery"). For small values of 7, this 
random code does not do as well as the other codes dis- 
cussed so far for the amplitude damping channel, but it 
still does significantly better than the case without error 
correction. Furthermore, for 7 > 0.35, our randomly gen- 
erated code actually outperforms all the other codes. For 
comparison, we have also plotted the worst-case fidelity 
for this randomly generated code in the absence of the 
transpose channel recovery, i.e., with the identity channel 
as the recovery map (line marked "random 4-qubit code, 
Id recovery"). One should keep in mind the ease with 
which the performance of the randomly generated code 
was achieved, due to the fact that the transpose channel 
is a near-optimal recovery map for any code. 



0.95 



fc 0.9 



0.8 




-[[5,1,3]] code 
-random 3-qul 

, random 2-qubit code, R p recovery 

■ random 4-qubit code, Id recovery 
-random 4-qubit code, R recovery 



0.05 0.1 



0.15 0.2 0.25 0.3 0.35 
damping parameter (y) 



0.4 0.45 0.5 



FIG. 2: Randomly generated two-, three-, and four-qubit 
codes using the transpose channel as the recovery map. For 
comparison, we have also plotted the worst-case fidelity for 
the [[5,1,3]] code, and that of the randomly generated four- 
qubit code with no recovery (i.e., identity channel as recov- 
ery). 

Finally, we also consider the possibility of construct- 
ing two-qubit and three-qubit codes for the amplitude 
damping channel. Because the transpose channel is near- 
optimal for any code, it can be used a good recovery map 
for the codes we generate, thus eliminating the need to 
search for a good recovery for every randomly selected 
code. The worst-case fidelity for the best codes we found 
are plotted in Fig. [2] For comparison, we also plot 
the worst-case fidelities for the randomly generated four- 
qubit code mentioned in the previous paragraph, with 
the transpose channel and the identity channel as recov- 



ery maps. The corresponding graphs for the two- and 
three-qubit codes with identity channel as recovery are 
close to that of the four-qubit code. From the figure, 
we see that while the worst-case fidelity decreases as the 
number of physical qubits decreases, the two- and three- 
qubit codes in fact do not perform too badly compared 
to the four-qubit code or the [[5,1,3]] code. Such codes 
may be of relevance whenever the desire to lower resource 
requirements trumps the need for the best possible worst- 
case fidelity. 



VII. CONCLUSIONS AND OPEN PROBLEMS 

In this work, we demonstrated the crucial role the 
transpose channel plays in perfect QEC, and used it to 
formulate a simple approach to characterizing and find- 
ing AQEC codes. Compared to previous work based 
on numerically-generated recovery maps specific to the 
noise channel in question, the universal and analytically 
simple form of our transpose channel makes it particu- 
larly useful towards developing a better understanding 
of AQEC. While not being the optimal recovery in the 
case of AQEC codes, the near-optimality of the trans- 
pose channel provides a simple algorithm for identifying 
codes that satisfy some maximum fidelity loss require- 
ments, without having to perform a difficult optimiza- 
tion over all recovery maps for every possible encoding. 
Furthermore, our approach, founded upon the worst-case 
fidelity rather than an averaged measure of fidelity, pro- 
vides the often desirable guarantee that the code found is 
able to protect all information that can be stored in the 
code with some minimum fidelity. We have also shown 
that the case of qubit codes is particularly easy to handle, 
and our method of computing the worst-case fidelity for 
a CPTP qubit map might be useful in contexts beyond 
our present discussion. 

There are many interesting related open problems. An 
immediate question is whether the gap present in our 
AQEC conditions between the necessary and sufficient 
conditions (arising from the inequality in Corollary [4]) 
can be reduced, either by improving the bound in The- 
orem [31 or by using a different recovery map that might 
perform better than the transpose channel. It would be 
very interesting if a simple and universal recovery map 
could be found, for which the dimension of the code does 
not appear in the worst-case fidelity. There is also the 
question of whether it might be possible to extend our 
efficient method of computing the worst-case fidelity for 
qubit codes to higher dimensional codes and more general 
channels. Finally, we expect that the transpose channel 
can also be used to study approximate codes more gen- 
eral than subspace codes, like for example, OQEC codes 

Another important problem is to figure out whether 
the transpose channel can be easily implemented us- 
ing measurements and gates. In the case of perfect 
QEC, the transpose channel (or cquivalently 7?.p e rf) can 
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be implemented simply using syndrome measurements 
and conditional gates (see for example, [22|). In or- 
der for AQEC codes to be useful for computational or 
communication tasks, it must be possible to implement 
the recovery operation using physical operations that are 
not overly complicated or demanding in resources. This 
is in fact another advantage of our analytical approach 
over numerically-constructed recovery maps for which no 
practical implementation structure may be apparent (al- 
though, see Q). 

AQEC provides a new and mostly unexplored arena of 
possibilities for the design of codes to protect informa- 
tion from noise for use in quantum information process- 



ing tasks. Our work provides an analytical characteriza- 
tion of AQEC and further analytical understanding will 
undoubtedly prove invaluable towards unlocking the full 
potential of AQEC. 
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